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Abstract 

A bigravity theory where the normal and dark matter components are cou- 
pled to separate metric fields linked to each-other with small non-derivative 
terms allows the Yukawa-like modification of the gravitational potential at large 
distances. This opens new prospects for the dark matter candidates. Namely, 
instead of being cold and collisionless, dark matter can be collisional and dis- 
sipative, as it occurs in the case of mirror matter that presumably does not 
form extended halos but is clumped similarly to visible matter. We show that 
the flattening of the galactic rotational curves can be explained if the typical 
scale of the Yukawa-like potential is about few tens of kpc and if the mass ratio 
between dark matter and visible matter in galaxies is about 10. 
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1 Introduction 



Cosmological observations show that the Universe is nearly flat with an energy density 
very close to the critical one: Qtot = + — 1- The dark energy component 

= 0.76 ±0.02 dominates today's universe while radiation is negligible and what is 
left is matter (ordinary+dark) with Qm = 0.24±0.02. In particular Qm = ^b + ^dm, 
where fls = 0.042 ± 0.005 and ^dm = 0.20 ± 0.02 are, respectively, the baryonic and 
dark matter components |T]. 

The standard cosmological model is in agreement with experimental results on 
Cosmic Microwave Background (CMB) and Large Scale Structure (LSS) that strongly 
support the presence of the dark matter component. Many DM candidates have been 
proposed in the literature axions {m ~ 10~^ eV), neutralinos (m ~ TeV), wimpzillas 
[m ~ 10^'* GeV) etc. However, the question why the fractions of the baryonic and 
dark components are so close, Qdm/^b ~ 5, remains unsolved. 

The most direct evidence for DM comes from galaxy rotational curves and cluster 
dynamics. In the Newtonian picture, far from the galactic center the velocity of 
rotating matter is expected to be of the form v{r) oc 1/y/r (known as Keplerian 
fall-off). This is due to the fact that the gravitational potential outside the galaxy 
inner core (bulge) is ~ GM/r, where M is total mass in the the bulge. Instead, 
one observes that outside the bulge v is approximately constant. The standard way 
to explain this anomalous behavior is to introduce non-baryonic dark matter (as in 
the standard CDM paradigm) distributed in an extended spherically symmetric halo 
around the galaxy, as a typical profile one can take an "isothermal" mass distribution 
p(r) oc (1 + (r/a)^)~^, where the scale a depends on the galaxj0. Since gravity 
is universal (no distinction between visible and dark matter), any point-like source 
(dark or ordinary) will generate the same potential 

0(r) = ^(Mi + M2) (1) 

where Mi and M2 is, respectively, the the mass of the visible and dark component. 
In the CDM scenario it is also assumed that the DM particles are coUisionless. As a 
consequence of this assumption, N-body numerical simulations tend to predict that 
DM clustering form singular profiles (cusps) p(r) oc l/r"", with a = l-i-1.5, [SllHllllE] 
quite different from the isothermal profiles. For most galaxies these cusp profiles can 
hardly reproduce the observed rotational curves [7], in particular, they provide for a 
steep growth of rotational curves at small distance. 

An alternative scenario is the mirror matter [S], according with our Universe is 
made of two similar gauge sectors. In other words, parallel to our sector of the ordi- 
nary particles (0-) and interactions described by the Standard Model, there exists a 

^The integral of the latter p{r) does not converge, therefore we need a cut-off radius L quite 
bigger than the visible galactic size. 
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hidden sector (M-) that is an exact duphcate of the ordinary sector with particles and 
interactions having exactly the same characteristics; the two sectors are connected by 
the common gravity (see for review [0]). Therefore, if the mirror sector exists, then the 
Universe should contain along with ordinary particles (electrons, nucleons, photons, 
etc.) also their mirror partners with exactly the same masses (mirror electrons, mir- 
ror nucleons, mirror photons, etc.). Mirror matter is invisible to us because it does 
interact with ordinary photons and naturally constitutes a dark matter candidate. 
One should stress that the fact that O- and M-sectors have the same micro-physics, 
does not imply that their cosmological evolutions should be the same. Indeed, if 
mirror particles had the same temperature in the early universe as ordinary ones, 
this would be conflict with Big Bang Nucleosynthesis (BBN). The BBN limit on the 
effective number of extra neutrinos implies that the temperature of the mirror sector 
T' must be at least about twice smaller than the temperature T of the ordinary sector 
allowing mirror baryons to be a viable candidate for dark matter. In particular, the 
mirror dark matter scenario would give the same pattern of LSS and CMB as the 
standard CDM if T'/T < 0.2 [lOJ. 

In addition, the baryon asymmetry of the Universe can be generated via out-of- 
equilibrium B — L and CP violating processes between ordinary and mirror particles 
whose mechanism could explain the intriguing puzzle of the correspondence be- 
tween the visible and dark matter fractions in the Universe, naturally predicting the 
ratio VLdm/^b ~ 1 ^ 10 [I2]. 

However, in contrast to the collisionless CDM, mirror baryons obviously constitute 
coUisional and dissipative dark matter. Therefore, one expects that mirror matter 
undergoes a dissipative collapse and thus DM (mirror) in galaxies has a distribution 
similar to the the visible matter instead of producing extended quasi-spherical CDM 
halos. Thus, for the DM sector as well for the luminous sector, an exponential profile 
p(r) oc e"''''^'" is expected. As a result, the distribution of the dark matter is more 
compact in the center of the galaxy and is not extended as in the CDM halos. If 
gravity is universal in the two sectors, the mirror dark matter hypothesis gets into 
difficulties in explaining the fiat rotational curves of galaxies. However, if each sector 
has its own gravitational field and there is a non-derivative mixing term [^, then 
gravity gets modified at large distance in way that rotational curves can be explained 
in terms of mirror DM. A detailed discussion of how gravity is modified is given 
in separate paper [131 El US]; for what concerns us the bottom line is that the 
non-derivative interaction between the two gravitational fields leads to a modified 
Newtonian potential. Namely, an ordinary test mass (matter of type 1) at distance 
r from a source of matter ordinary + dark mater (matter of type 2) (Mi and M2), 

^The particle mixing phenomena between ordinary and mirror sectors were discussed in the 
Hterature for photons J£J, neutrinos p^, neutrons [TH], etc., as well as possible common gauge 
interactions between two sectors [H]. 
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instead of ([T]) is the subject to the following potential 



0(r) = ^ (Ml + M2) + %^ (Ml - M2) , (2) 
2r 2r 

where G is the Newton constant and is the graviton inverse mass. Notice that 
at small distance r <^ gravity is not modified: an ordinary particle feels only the 
influence of luminous ordinary matter Mi through the standard Newton potential: 

Ml 

ct)(r < Tg) ^ (3) 

At large distance r ^ rg the interaction with dark matter kicks in and effectively the 
particle "feels" the sum of all matter Mi + M2 with a halved Newton constant G: 

, /G'\Mi + M2 
<P{r » Tg) ^ J \ (4) 

This is the key result that enables us to reproduce the observed rotational curves of 
galaxies. 



2 Galaxy Rotational Curves 

A rotational curve describes the velocity of gravitionally bounded galactic objects 
(typically stars and interstellar gases) as a function of the distance r from the center. 
Most of matter a galaxy (e.g. a disk galaxy) is concentrated in the internal region, 
called bulge. In the simplest schematization, a galaxy is modelled as an inner spher- 
ically symmetric constant density region; the picture can be refined taking a matter 
density profile that follows approximately the luminosity profile which decreases ex- 
ponentially with r moving from the center to the outer region: 

p(r)=Poe-^/-, Po = £-s, (5) 

where r = x"^ + y"^ + z"^ (the galactic plane lies in the xy plane), m is the the galaxy 
total mass and po is the bulge size. The total visible mass of a galaxy can be estimated 
from the the bulge size po and the velocity (wq) in the region where the plateau starts 
as 

In a different cylindrical symmetry modellization the galaxy is taken to be two- 
dimensional with a disc surface density 

a(r.) = aoe-/-, ao = £^2^ (7) 
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Here Vc = \/ + y'^. A fully realistic modellization of the matter distribution is an 
intermediate between the spherically and cylindrical symmetric cases. For instance, 
due to the rotational symmetry along the z axis, the exponential fall-off will depend 
not only by the distance from the center but also on z 



Pir) = Pfoe -of{z/zo) 



(8) 



where tq and zq are of order of the bulge size (we will refer to it as "f- modulated"). 
When / = 1 we have spherical symmetry. In general, a good approximation is 
f{z/zo) = sedi^iz/zo). 

Consider now the radial velocity distribution t>(r) for matter (gas, stars, etc.) 
rotating around the galactic center. According the Newtonian theory, v{r) is expected 
to fall-off as 1/ y/r at large distances. Observations show a totally different behavior: 
asymptotically, there is a plateau (flattening). In the standard paradigm to reconcile 
theory and experiments a new kind of "dark" matter is introduced providing galaxies 
with an invisible spherically symmetric halo. The rotational velocity v"^ has then two 
contributions: one from the visible matter distributed according to ([H]) and the second 
from the dark component. 

(9) 



The distribution profiles strongly depend on the nature of dark matter, a common 
parametrization is 
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(10) 



where is a characteristic length scale. The values for a,/?, 7 are shown in the 
table (IT=Isothermal, NWF=Navarro-Frenk- White and M=Moore). By using the 
IT model one can nicely fit many galaxy rotational curves, even those dominated by 
the dark matter as the dwarf and low surface brightness (LSB) galaxies. The trouble 
with IT model is that, when N-body simulations are used to predict DM profiles, 
singular distributions are produced, containing a cusp when r ^ 0. For instance in 
the NFW model [21 [3] p ~ 1/r, while in the M model ^ i5j p ~ l/r^-^ (see also [6]). 
These profiles are shown in Figures. ([T]) and ([2]). Figures ^ and (jl]) compare the 
visible matter contribution to the total velocity for profiles with different symmetry. 
In the cylindrical case the bump at low distance from the center is enhanced. Finally, 
in figures ( fT5|) and ([5]) the visible and dark contribution in the IT model both for 
spherical and cylindrical symmetric profile. Let us now discuss how the rotational 
curves are reproduced using the potential In order to grasp basic the idea, let us 
imagine a galactic object moving along a circular orbit of radius r around the center; 



5 




Figure 1: Density profiles for: IT (solid), 
NFW (dashed) and M (dotted) in arbitrary 
unit. 
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Figure 3: Visible matter contribution to 
the rotational curves: spherical (solid) and 
cylindrical approxiamtion (dashed), with 
ro = 3 kpc and M = IO^^Mq. 



Figure 2: Rotational curves (Spheri- 
cal approximation) for: IT (solid), NFW 
(dashed) and M (dotted) in arbitrary unit. 
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Figure 4: Visible matter contribution to 
the rotational curves: sphericaal (solid) and 
f-modulated approximation (dashed), with 
ro = 3 kpc and M = W^^Mq. 



its velocity is determined by equating the centrifugal force with the radial component 
of gravitational pull g{r) derived from the potential ([2]) 



f^(r) 



g{r) = G 



Ml + M2 _^ Ml - M2 



IH I e '■f 

r, 



11) 



2r2 2r2 

Notice that if Mi = M2 we are back to the Newtonian theory; on the other hand 
when Ml 7^ M2 and r ^ we get 

'G\ m(Mi + M2 



9m[r) 



(12) 



effectively we have a Newtonian force with the Newton constant G replaced by G/2. 
When r <^ Tg a ordinary test mass m will feel a force only from visible matter sources 

mMi 



G- 



(13) 



Summarizing, when r << rg, the influence of dark matter on visible matter is negli- 
gible and the behavior is Newtonian. In the opposite limit: r >> rg, the force is also 
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Figure 5: Different contributions to the 
rotational curves in IT dark matter halo 
(cylindrical approximation): (dashed), in- 
visible (dotted) and total (solid), with vq = 
8 kpc, Vy = 3 kpc and Mdm = 5My. 
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Figure 6: Different contributions to 
the rotational curves IT dark matter 
halo (spherical approximation): visible 
(dashed), invisible (dotted) and total 
(solid), with ro = 8 kpc, r„ = 3 kpc and 

Mdm = 5My. 



Newtonian, though a test particle feels the total mass Mi + M2 and the effective New- 
ton constant is G/2. In the crossover region r Tg there is a sizeable deviation from 
Newtonian theory due to the presence of type 2 matter, resulting in a enhancement 
of f ^ as shown in Figs. ([7]) and ([H]) for point-like distributions. 
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Figure 7: Rotational curves v{r) for point- 
like distribution: Newtonian (solid) and 
modified potential for M2 = lOMi and 
rg = 15 kpc (dashed). The enhancement in 
the shape for modified gravity can explain 
the observed flattening. 
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Figure 8: Normalized point-like rotational 
curves v{r)r^/'^ (M2 = lOMi and rg = 15 
kpc). The flat behaviour at large scale is in 
fact a slow fall-off as r""*^/^. 



For extended objects one has the following integral expression for the acceleration 

;%)^(r-,)(l + ^)e-^l. (1 

2|r — r l'^ \ "^g J J 
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where pi,2(r) are the density distribution for the matter of type 1 and 2, respectively, 
r is the position of the test mass. For the benefit of reader we have collected in a 
table the relevant expressions for different profiles discussed. 
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Here, Mi and M2 are the total mass of the visible and dark component of the galaxy 
respectively; ri^2 and 2:1,2 are their bulge sizes and p/1,/2 = (/ d'^r' fi^2iz/ zo)e''^^^'^)~^ . 
For both the spherical and cylindrical symmetry case, '& is the angle between r and 
r'. By symmetry, the only non- vanishing component g{r) of g is along r. Therefore, 
the circular rotational velocity, v{r), is given by 

f (r) = \/rg{r) (15) 

Though for the potential ^ Gauss theorem does not apply, the deviation is not 
numerically relevant as one can see from Fig.©. In the model presented it is natural 
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Figure 9: Comparison between rotational 
curve (spherical) computed with numeri- 
cally method (solid) and using the Gauss 
theorem (dashed). The result is similar 
(M2 = lOMi, ri = 3 kpc, r2 = 5.4 kpc, 
= 15 kpc). 
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Figure 10: Rotational curves with modi- 
fied potential: comparison between spher- 
ical (dashed) and f-modulated (solid) ap- 
proximation {vg = 20 kpc, M2 = lOMi, 
ri = 3 kpc and r2 = 5.4 kpc). 



to consider that matter of type 1 and type 2 distributed in the same way. Indeed, each 
type of matter have the same non-gravitational interactions and the sectors 1 and 2 
couple only gravitationally. The very same formation mechanism for visible galaxies 
(type 1) will also form the invisible (type 2) ones giving the same matter distributions 
up to a rescale of relevant parameters . The scale R2 is set simply by imposing that 
p2 is related to pi by a dilatation factor a fixed by the invisible to visible mass 
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ratio (3 = M2/M1. Then it follows that R2 = aRi, with a = (Ms/Mi)^/^ = (3^/^. 
For instance, if i?i = 3 kpc and (3 = 10, then R2 ~ 5.4 kpc. For the cylindrical 
approximation, a similar discussion tell us that a = \/]3. 

As a result, matter of type 2 plays the role of the dark matter with an additional 
welcome benefit: there is no clash between the actual halo shape (i.e. that fits the 
observed rotational curves) and the dark matter density profile produced by numerical 
simulations. The invisible, dark, matter of type 2 has a distribution similar to visible 
one and the observed rotational curves is reproduced by the non- Newtonian character 
of the gravitational potential at intermediate distances. 

A large analysis varying the main parameter of the theory are shown in Figs. (fTOl) . 
(El), (III), (III, (El), (El) and (IIHl^ Both variations of M2/M1 and can determine 
the flattening of the rotational curve for a large range of distances. Another relevant 
term is the sizes of the bulges of both matter distribution. It is interesting to note 
(see. Figs. (fT^ and (|TU|) ) that the bigravity rotational curves can approximately 
reconstruct a shape similar to the shape obtained with the classical IT model. Notice 
that our model requires dark matter two times more abundant than the matter needed 
in the CMD model. This amount is the same even for cluster dynamics and cosmology 
because, at large distance, gravity is Newtonian and reproduces the same results of 
the standard CDM with Mdm ^ Mb- 

The value of the parameter Vg stays in the range 10 -r- 20 kpc, that is the typical 
size of a galaxy (M ~ 10^^ M©). 
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Figure 11: Different behaviors of ro- 
tational curves (Spherical approximation) 
varying M2/M1 mass ratio: 1 (dotted), 5 
(dashed) and 10 (solid) (r^ = 20 kpc, ri = 3 
kpc and r2 = 5.4 kpc). Notice that, for 
Ml = M2, we obtain once again the New- 
tonian shape. 
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Figure 12: Different behaviors of ro- 
tational curves (Spherical approximation) 
varying the r^: 10 kpc (dotted), 20 (dashed) 
and 40 (solid) (M2/M1 = 10 kpc and n = 
10 kpc). 



Finally, we show a comparison with the experimental data relative to the NGC 
2403 galaxy (see Figs. ( |T9l) and (l20l)). Observe that in this disk galaxy the cylindrical 
shape approximation reproduce with high accuracy the profile observed, while the 
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Figure 13: Comparison between spherical Figure 14: Rotational curves with mod- 



(dash) and cylindrical (solid) approxima- 
tion in modified potential scenario. Here, 
n = 3 kpc, r2 = 5.4 kpc, Mi = 1, 
Ma /Ml = 8,rg = 20 kpc. 




Figure 15: Comparison between rota- 
tional curves (spherical approximation): 
modified potential (dashed) (M2 = lOMi, 
ri = 3 kpc, r2 = 5.4 kpc and rg = 20 kpc) 
and IT dark matter halo (solid) (r^ = 8 kpc, 
r„ = 3 kpc and Mdm = 5Mu); Keplerian fall 
(dotted). 



ified potential: comparison between cylin- 
drical (dashed) and f-modulated approxi- 
mation (sohd) {rg = 20 kpc, M2 = lOMi, 
ri = 3 kpc and r2 = 5.4 kpc). 
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Figure 16: Comparison between rota- 
tional curves (cylindrical approximation): 
modified potential (dashed) (M2 = lOMi, 
ri = 3 kpc, r2 = 6.5 kpc and rg = 20 kpc) 
and in IT dark matter halo (solid) (r<j = 8 
kpc, r^ = 3 kpc and Mdm = 5M„). 



spherical approximation underestimates the values of velocity at low distances. 

3 Discussion and Conclusion 

Since DM is more abundant will start collapsing before ordinary matter and the 
latter will fall in the potential well of DM. The dark mirror matter is coUisional 
and, differently from CDM, undergoes a dissipative collapse creating structures (e.g. 
galaxies, etc.) similar to the visible ones, i.e., with the same mass distribution. Since 
the gravitational force between sector 1 and 2 is effective when r > ~ 10 -i- 20 
kpc, after an initial stage of mutual interaction, as the collapse proceeds, the inner 
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Figure 17: Different beliaviors of rota- 
tional curves (Cylindrical approximation) 
varying M2/M1 mass ratio: 1 (dotted), 8 
(dashed) and 15 (solid) (r^ = 20 kpc, ri = 3 
kpc and r2 = 6.4 kpc). Notice that, for 
Ml = M2, we obtain once again the New- 
tonian shape. 
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Figure 19: Fit (cylindrical approxima- 
tion) NGC 2403: n = 3 kpc, rs = 6.5 kpc. 
Ml = 0.26 1O^^M0 , M2/M1 = 10, Tg = 10 
kpc. 
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Figure 18: Different behaviors of rota- 
tional curves (Cylindrical approximation) 
varying the r^: 10 kpc (dotted), 15 (dashed) 
and 20 (sohd) (M1/M2 = 10 kpc and ri = 3 
kpc). 
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Figure 20: Fit (spherical approximation) 
NGC 2403: n = 2 kpc, rs = 10 kpc. Mi = 
0.26 IO^^Mq , M2/M1 = 10, rg = 10 kpc 



structures of DM and ordinary matter are formed independently. On the contrary 
one can argue that for the outer region of galaxies, the interaction between ordinary 
and dark matter is important in the angular momentum exchange giving rising to a 
common galactic plane. 

Besides the flattening of the rotational curves, a good DM matter candidate should 
also reproduce some properties that have observed in galaxies of different size and 
type PE]. 

Typically, to fit the DM matter profile a constant inner core is required. In our 
model, dark and visible matter follow the same distribution and an inner core region 
with an exponential profile p{r) oc e~^^^° {r^ is the DM bulge size) naturally exists 
like for the visible galaxy. Moreover, in the standard paradigm, to fit the data it is 
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assumed that DM is present only in a region with r < L, and the cut-off radiuqj L is 
given by the following phenomenological law L = 1OO(Mdm/1O^^M0)^/^ kpc. 

The latter can also be obtained in our model if make the rather natural assumption 
that L oc r£,. Indeed, we have an exponential DM distribution, with total dark mass 
given by Md = pDSnr j^ as a result L oc M^^ . 

Finally, observations show that in a galaxy the DM amount is inversely propor- 
tional to the visible mass. This property is not surprising in our model where the 
graviton Compton length is of order 10 ^ 20 kpc: for scales r > Vg Gravity is 
Newtonian and in a galaxy with an optical radius greater than 30 -^ 40 kpc the rota- 
tional curve shows that resembles the Keplerian fall; this is interpreted in standard 
paradigm as the observation that bigger galaxies have less DM. 

Intriguingly, recent astrophysical observations indicates for dark matter a peculiar 
behaviour in galaxy cluster collisions. In the Bullet cluster ^20J, the dark matter 
component of the cluster shows a coUisionless behaviour; instead, in the Abell 520 
cluster [21] there is some evidence for an important contribution from self-interaction. 
In mirror dark matter models those observations are not surprising [22] . 

In conclusion if the Universe is made of two separate sectors, one visible and one 
hidden, each of them can have its own gravity and the two metrics can interact leading 
to a large distance modification of the gravitational force. At small distances, where 
the gravitational interaction between the different sectors is effectively shut off, the 
behaviour is Newtonian. The mirror matter as dark matter, together with Lorentz 
breaking can explain the flattening of galaxy rotational curves. This model assumes 
dark matter and modification of gravity. According to this proposal dark and visible 
matter have mass distribution avoiding the tension between observed rotational curves 
and extended CMD halo distributions predicted by N-body numerical simulations. 

It is worth to point out that, at large cosmological distance both sectors interact 
each-other with effective Newton constant G/2. It follows that the observed Hubble 
constant implies that the total energy density of the Universe is twice bigger than 
in the standard cosmology in which the Newton constant at the same cosmological 
distance remains the canonical G. In this scenario, instead of pcr = 3Hq/8ttG implied 
for the standard cosmology, we must have pcr = 3Hq/AttG. 

The mirror matter model with modified potential does not exclude the possibility 
of direct (non-gravitational) interaction between the normal and dark matter com- 
ponents, for example via photon kinetic mixing [13]. This also makes possible the 
direct detection via such interaction [23] and can explain the DAMA results [21]. In 
our model this possibility becomes higher because of the twice bigger DM amount 
in galaxy and the related bigger velocity of dark particles. Let us remark also that 
our scenario can have interesting implications also for the neutrino-mirror neutrino 
(active-sterile) oscillations [IT] and can make more attractive the case of the exact 

■^For IT models cut-off is needed to have a finite the total DM mass. 
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mirror parity when the ordinary and mirror neutrinos are degenerate in mass and 
maximally mixed. But their oscillations inside the galaxy should suffer a MSW-lilce 
suppression because of different gravitational potential felt by active and sterile neu- 
trino species. Far outside the galaxy where the gravity becomes universal oscillations 
instead can go in full strength. 
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